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Abstract 

The purpose of this paper is to present a comparison between the modified nonlinear Schrodinger 
(MNLS) equation and the focusing and defocusing variants of the (unmodified) nonlinear Schrodinger 
(NLS) equation in the semiclassical limit. We describe aspects of the limiting dynamics and discuss 
how the nature of the dynamics is evident theoretically through inverse-scattering and noncommutative 
steepest descent methods. The main message is that, depending on initial data, the MNLS equation can 
behave either like the defocusing NLS equation, like the focusing NLS equation (in both cases the analogy 
is asymptotically accurate in the semiclassical limit when the NLS equation is posed with appropriately 
modified initial data), or like an interesting mixture of the two. In the latter case, we identify a feature 
of the dynamics analogous to a sonic line in gas dynamics, a free boundary separating subsonic flow from 
supersonic flow. 



1 Introduction 

This paper is concerned with the semiclassically scaled modified nonlinear Schrodinger (MNLS) equation: 

le S + r 2 S + l ^ + m£ cl (l ^ )=0 ' a ' £>0 ' W 

which can be thought of as a perturbation (with perturbation parameter a) of the nonlinear Schrodinger 
(NLS) equation: 

taken in the so-called focusing case of n = +1. (For k = — 1 one has the so-called defocusing case.) Our 
goal is to describe connections between the MNLS equation on the one hand and the NLS equation (perhaps 
surprisingly, both focusing and defocusing types are related to MNLS) on the other, connections that become 
especially clear in the semiclassical limit e<l. 

*The work of the first two authors was partially supported by the National Science Foundation under grant DMS-0807653. 
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1.1 Formal considerations. 



Erwin Madelung was probably the first person to propose thinking of Schrodinger-type equations for a 
complex field <f> = <p e (x,t) as describing the motion of a kind of fluid [TU]. He introduced the real- valued 
fluid dynamical variables p e (density) and u e (velocity) given in terms of 4> e by 

p e {x,t) := \cj) t (x,t)\ 2 and u e (x, t) := e3 | ^- \og{(j) e {x, i))| . (3) 

Here \og((f) t (x : t)) is defined by continuation to be a smooth function of x. A key observation is that if the 
complex- valued field <p t {x,t) initially has the form of an oscillatory wavepacket: 

e (x,O) =A(x)e lS{x)/t (4) 

where A and S are real-valued functions independent of e, then the corresponding initial values of the 
Madelung variables are independent of e as well: 

Po{x) := Pe{x, 0) = A{x) 2 and uo(x) :— u e (x,0) — S'(x). (5) 

Without approximation, the NLS equation |2]) can be written in terms of the Madelung variables as: 

c>Pe . d , , n du e d fl 2 \ 1 2 dF[p e ] 

-dt + dx {p ^= and -w + dx [2 u '- Kp ' rr^r (6) 



where F[p] denotes the expression 

P\n\ ■ = 

2p dx 2 \2p dx J 

With e-independent initial data of the form ^ taken in it seems reasonable to neglect the dispersive 
term ^e 2 d x F[p e ] to obtain limiting e-independent equations that one might expect to govern the leading 
terms (p(x,t),u(x,t)) of the e-dependent fields (p e (x,t),u £ (x, t)) when f«l. The e = truncated system 
of equations 

dp d , . , flu 9 /l , \ 

^ + -M = and g t +Q x {-^-^)=0, (8) 

which we call the dispersionless NLS system, is really a classical Euler system of compressible fluid (gas) 
dynamics. The term — up in the flux for u represents the pressure as a function of density. If k = — 1 
(the defocusing case of the NLS equation) then the pressure is an increasing function of the gas density 
and the physical dynamics is that of the gas "trying to get out of its own way" ; localized concentrations in 
density tend to decay with time. On the other hand, if k — +1 (the focusing case of the NLS equation) 
then the pressure decreases with increasing density and localized concentrations in density are dynamically 
enhanced. In physical problems with such pressure laws, like supercooled Van der Waals gases, initially 
small but nonuniform perturbations of a uniform density lead quickly to spontaneous condensation of the 
gas into liquid droplets. Mathematically, the equations Q comprise a quasilinear system of hyperbolic type 
for k = — 1, and the Cauchy initial- value problem for such systems is locally well-posed in suitable function 
spaces. On the other hand when k — +1 we have instead a system of elliptic type, and the Cauchy problem 
is ill-posed. 

The ill-posedness of the Cauchy initial-value problem for the e = system ^ in the focusing (k = +1) 
case can be regarded as a manifestation of the well-known modulational instability of the focusing NLS 
equation. This instability becomes increasingly severe for e <C 1 and is asymptotically catastrophic in the 
semiclassical limit e 4- 0. The modulational instability can lead to unwanted physical effects in systems 
modeled by the focusing NLS equation. For example, the Gordon-Haus jitter effect of pulses in optical fiber 
telecommunications systems is caused at least in part by the modulational instability. In such situations, 
physical effects originally neglected to derive a simple model equation (the focusing NLS equation) might 
play a role in reducing the influence of the instability. One effect that can be included perturbatively is 
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nonlinear dispersion, and including such a term results in the MNLS equation ([lj in place of the focusing 
NLS equation. Obviously if a — the MNLS equation becomes the focusing NLS equation, but for a > 
the dynamics can be different, perhaps usefully so. 

In terms of the Madelung fluid-dynamical variables the MNLS equation ([I]) becomes 

dp, d f 3 ,\ „ , du t d [1 2 \ 1 2 dF[ Pe ] 

where again F[p] is given by Q. The formal limit system (set e = 0) is in this case the dispersionless MNLS 
system 

d P , 9 , 3 _2\ , 9U , 9 /In 



This quasilinear system is obviously an a-perturbation of the corresponding elliptic-type dispersionless NLS 



system <[8j) with n = +1. Unlike the a = case, the system ( 10 1 is of mixed type, being hyperbolic for Q > 
and elliptic for Q < 0, where 

Q := a 2 p + au — 1. (11) 

This simple calculation suggests that the modulational instability of the focusing NLS equation (a = 0) 
can be completely suppressed by nonlinear dispersion (a > 0) for initial data satisfying the condition 
a 2 po(x) + auo(x) — 1 > 0. Quasilinear systems of mixed type occur frequently in problems of transsonic 
gas dynamics. In regions of spacetime where the gas velocity is subsonic (respectively supersonic), the 
equations are of elliptic (respectively hyperbolic) type. Throughout this paper we will use this terminology 
in the context of the mixed-type system (|10[) also. To make this analogy more obvious, we can define a 



density-dependent sound speed for the dispersionless MNLS system ( 10 1 by 



c(x,t) := --ap(x,t), (12) 
a 

and then recalling that u(x, t) has the interpretation of fluid velocity, the condition of hyperbolicity for the 



system (10 1 becomes the supersonic condition u(x,t) > c(x,t) while that of ellipticity of the system (10) 



becomes the subsonic condition u(x,t) < c(x,t). 



The fact that the dispersionless MNLS system ( 10 ) can be either hyperbolic like the dispersionless de- 



focusing NLS system ((|8| with k = —1) or elliptic like the dispersionless focusing NLS system ((|8j) with 
K = +1) depending on the sign of Q establishes a qualitative connection between these systems in the semi- 
classical limit. However, the connection is far more concrete as was shown in 2J. Indeed, let us suppose 
that p(x,t) and u(x,t) are functions satisfying the dispersionless MNLS system ( [To| and the (supersonic) 
hyperbolicity condition u(x,t) > c(x,t). Then if new fields p(x,t) and u(x,t) are defined in terms of p(x,t) 
and u(x, t) by the algebraic mapping 

p := pQ and u := u + lap for Q > 0, (13) 

it is easy to check that the pair (p, u) constitutes (after dropping the "hats" ) a solution of the dispersionless 
defocusing NLS system with k = —1). Similarly, if p(x,t) and u(x,t) solve the dispersionless MNLS 



system (10 1 and satisfy the (subsonic) ellipticity condition u(x, t) < c(x,t), and if new fields p(x,t) and 



u(x, t) are instead defined by 

p := — pQ and u := u + 2ap for Q < 0, (14) 
then (p,u) constitutes a solution of the dispersionless focusing NLS system ( dsb with k = +1). Note that 



both mappings (13) and (14) preserve the positivity of the Madelung fluid density p. This simple calculation 
reveals the remarkable fact that, in the semiclassical limit, the dynamics of the MNLS equation can be 
exactly mapped onto either the dynamics of the defocusing NLS equation or those of the focusing NLS 
equation depending only on the sign of Q. 
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1.2 Inverse spectral theory. 



One of our aims in writing this paper is to show how the formal observations made above can be placed 
on completely rigorous mathematical footing, and also how the interesting transsonic case where Q is sign- 
indefinite can be analyzed with great precision in the semiclassical limit. The main tool here is the exact 
solution of the Cauchy initial-value problem for the MNLS equation (JlJ by means of the inverse spectral 
transform specially adapted to that equation. We also wish to make a comparison with the solution of the 
NLS equations of both focusing and defocusing types, which is carried out by a better-known inverse spectral 
transform. 



1.2.1 Inverse spectral transform for MNLS. 



For the reader's convenience, we give here a very brief description (adapted from the appendix of [5]) of 
the calculations involved with implementation of the inverse-spectral transform for the MNLS initial-value 
problem. Let <po = 4>{x) denote the complex-valued initial condition for the MNLS equation, a function that 
may depend parametrically on e > as in the wavepacket case Q and that is assumed for convenience to 
be of Schwartz class. Let k be a complex spectral parameter. Given <j) , the Jost matrices 3±(x;k, e) are 



(uniquely) defined for k G 



to satisfy the linear differential equation 
dJ± 



where the coefficient matrix L is given by 

L := Acr 3 + 2ik 

and where A is defined by 



dx 





t>o(x)* 



LJ. 



<f>o{x) 






'1 " 




-1 



A 



2i 



(15) 



(16) 



(17) 



as well as the asymptotic normalization conditions 3±(x;k,e)e Axcr 3/e j as x _^ ±oo. The Jost matrix 
solutions are unimodular for all x G R and k 2 € R, and in particular are fundamental matrix solutions of 



(151. Since the differential equation (15) can have only two linearly independent column vector solutions, 
the columns of 3+{x]k, e) are linear combinations of those of 3-(x;k,e) and hence there exists a 2 x 2 
unimodular scattering matrix S(fc, e) defined for k 2 G R by 



J+(x; fc, e) = J_ (x; k, e)S(k, e), k 2 G 



(18) 



In particular, the ratio 

Su(k,e) 
S22(k,e)' 

is called the reflection coefficient corresponding to the potential 4>q; r(-,e) should be viewed as a kind of 



k 2 e 



(19) 



transform of 4>q(-) derived from the behavior of solutions of ( 15 1 as k varies through the continuous spectrum. 
By comparing the Jost solutions for — k with those for k it can be shown that r is an odd function of k. 

There is in general also information about 4>a encoded in the discrete spectrum of the problem (15 1. It 
turns out that the second column of J+ and the first column of J_ can be analytically continued into the 
region 3{fc 2 } > 0. It follows that 5 l 22(fc, e) extends analytically to the region 9{/c 2 } > 0. If kj G C with 
Q{k 2 } > is a zero of S22{k, e), then it is easy to see that (x; kj, e) and j_(x; kj, e) are proportional, and 
hence there exists a nonzero constant 7^ G C such that 



j+ ' [x] kj , 



Tjj-^fcj^), 



5 22 (fc„e)=0, %{k 2 }>0 



(20) 



holds for all x G R. Such k = kj are eigenvalues or discrete spectrum of the differential equation ( 15 ) 



because it can be shown that the condition of proportionality (|20|) implies the existence of a solution of (|15|) 
for k = kj that is exponentially decaying both as x — > 



and also as 
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If (f>o is such that it turns out that S22{k, e) has a finite (for fixed e > 0) number of zeros in its domain 
of analyticity, all of which are simple zeros (so that S22(kj,e) ^ 0, where prime denotes differentiation 
with respect to k), and that S22(k) does not vanish at all for 5{fc 2 } = 0, then the corresponding solution 
4> = 4> c (x,t) of the MNLS equation with initial condition (f> a can be obtained from the scattering data 
consisting of the reflection coefficient r(k, e) for 3{fc 2 } = together with the pairs (kj,jj) corresponding to 
the discrete spectrum by the solution of a matrix Riemann-Hilbert problem. 

The Riemann-Hilbert problem can be formulated to take advantage of a certain symmetry of the differ- 
ential equation ( 15 1 under the reflection k n- — k. Thus we introduce the spectral variable z = k 2 and seek a 
2x2 matrix N(z) = N(^; x, t, e) satisfying the following conditions (N(z) := M^z 1 / 2 ), where M(/c) denotes 
the matrix solution of the Riemann-Hilbert problem as formulated in the appendix of [5] and where z 1 / 2 
denotes the principal branch of the square root). Let 



and for z € R set 



9(z; x, t) 



s(z,e) := r(e i7r/4 (-^) V V) 



2 














H) 


a 




or 





(21) 



x l 2 e) = 



(22) 



Finally let D C C denote the finite set D = {k\, . . . , k 2 Nl k\ 2 
the following conditions: 



r(z 1 / 2 ,e), 
r{i{-z) 1 '' 1 ,e), 

, fc^}. We seek a matrix N(z) satisfying 



z > 
z < 0. 



Analyticity: N(z) is an analytic function of z for z € C\ (lUfl), taking continuous boundary values 
N±(z) on the real line from C±. 



Jump condition: The boundary values taken by N(z) on R are related as follows: 

1 



and 



N+(z) = N_(z) 
N + (z) = i ff3 N_(z)r <T3 



s(z,e)*e~ 2ie< - z:x ^/ e 
1 

s{z,e)*e- 2l9 ^ x ^/ e 



-s(z,e)e 2l9( - z ' x ^/ t 
l + \s(z,e)\ 2 



l-\s(z,e)\ 2 



z > 0, 



z < 0. 



(23) 
(24) 



Point singularities: The matrix N(z) has simple poles at the points of D. If Zj = k 2 E D with 
3{zj} > 0, then 



_ c , e 2i9(zj\x,t)/e 





-2i9(z*;x,t)/e q 



(25) 



Res N(z) = lim N(z) 
Res N(z) = lim N(z) 

z—z* z—tz* 
J J 

where Cj := — Zkjjj / S^ikj, e). 

Normalization: The matrix N(z) is normalized to the identity at the origin in the complex plane: 
N+(0) = N_(0) =1. 

From the solution of this Riemann-Hilbert problem, one extracts the solution (f) e {x,i) of the Cauchy initial- 
value problem for the MNLS equation ([!]) with initial data </) e (x, 0) = <fio(x) from the formula 



a"™ N22\z;x,t, e) 



(26) 



(Note that N(z) = N + Niz' 1 / 2 + C^z" 1 ) as z — >• oo, where No and Ni are well-defined as functions of x 
and t, and moreover Nq is diagonal [5J page 990] with det(No) = 1.) 
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1.2.2 Inverse spectral transforms for focusing and defocusing NLS. 

As was first shown by Zakharov and Shabat |18j . the initial- value problem for the NLS equation ^ with 
Schwartz-class initial data 4>o can be solved by means of the scattering problem 



dJ 



ZS,K 



ZS,K t ZS,k 



dx 



where the coefficient matrix L ZS,K is given by 



• ZS,K 





-K(j) (xy 



(f> (x) 




(27) 



(28) 



where A € R is a spectral parameter. Here as before J ± ' 
boundary conditions that 3^' K (x; X,e)e iXa3x/e -> I as x 
matrix S ZS ' K (A, e) by the relation 



ZS,K 



= 3±' K (x;A,e) are Jost matrices satisfying the 
±oo, and once again one defines a scattering 



jf' K (x; A, e) = 3^' K (x; A, e)S zs ' K (A, e), A e 



,ZS 



and a reflection coefficient by 



.zs. 



K (A,e) 



A G 



(29) 



(30) 



This completes the description of the scattering data associated with the continuous spectrum of the direct 
scattering problem (27). 

If re = — 1 (the defocusing case), then as an eigenvalue problem with eigenvalue A, the scattering problem 
(27) is self-adjoint for </>o Schwartz class, and there can be neither complex eigenvalues nor embedded real 
ones. But if re = +1 (the focusing case), then selfadjointness is lost and there can be discrete spectrum. 



To find the associated scattering data, note that the second column of J ZS ' K and the first column of J 



,ZS,K 



extend analytically into the upper half-plane 3{A} > 0, and hence so does S 22 ' re (A, e). While for re = — 1 
the latter function is bounded away from zero in the closed upper half-plane, if re = +1 it may have zeros, 
which are the discrete eigenvalues of the problem p7| for re = +1. If in the latter case Aj G C+ is such a 



(31) 



zero, then the corresponding columns of the Jost matrices are proportional for A = Aj 



•zs,«(2), . 



^22 '"(Aj, e) 



0, 9{A 3 }>0, re = +l 



for some nonzero constant jf S - 

In the defocusing case, the solution of the Cauchy initial-value problem with Schwartz-class initial data 
00 can always be obtained from the reflection coefficient alone via the solution of a matrix Riemann-Hilbert 
problem. In the focusing case we must additionally assume that </>o is such that the number of zeros of 



£22 (A, e) in the closed upper half-plane is finite, and all are nonreal and simple, and we must also include 
the pairs (Xj, 7 ZS ) as part of the scattering data. Then, setting 



(A; x, t) := —Xx — X 2 t, 



(32) 



and letting D zs :— {Aj., . . . , Ajv, A*, . . . , A^} denote the eigenvalues and their complex conjugates (this 
is an empty set in the defocusing case), the Riemann-Hilbert problem is to find a matrix N ZS,K (A) = 
N ZS,K (A; x, t, e) satisfying the following conditions: 

Analyticity: N ZS,K (A) is an analytic function of A for A G C\ (KUD ZS ), taking continuous boundary 
values N^. ' K (A) from C±. 

Jump condition: The boundary values taken by N ZS ' K (A) on K are related as follows: 



N Z _ S ' K (A) =N^' K (A) 



ZS,k / 



„ZS,K 



-nr 



ZS.K 



(A,e)*e~ 2lf,ZS(A;:M)/e 



(A,e)e 



2i6 zb (\;x,t)/e 



l + re|r /fa ' K (A,e)| 2 



z G 



(33) 



G 



Point singularities: The matrix N ZS,K (A) has simple poles at the points of D zs (an empty set for 
k = -1). If Xj £ D zs with 3{Aj} > 0, then 



Res N (A) = lim N ,K (A) 

A — A i A — >Ai 



Res N ZS < K (A) = lim N ZS ' K (A) 

A=A* A-vA* 



_ c ZS e 2i8 zh (\ j ;x,t)/e 



c ZS* e -2i9 zs (\*;x,t)/e q 



(34) 



where cf := - 7 f /^'(A,, e). 



Normalization: The matrix N ZS ' K (A) is normalized to the identity at infinity in the complex plane: 
N ZS < K (A) -> I as A -> oo. 

From the solution of this Riemann-Hilbcrt problem, one extracts the solution <p e (x,t) of the Cauchy initial- 
value problem for the NLS equation ([2| with initial data <j> e (x, 0) = (f>o(x) from the formula 



(j) e (x,t) :=2i lim XNf 2 ' K (X;x,t,e). 

A— foo 



(35) 



1.3 Outline of the paper. 

In this paper we will explain how in the semiclassical limit the MNLS equation can be regarded as a 
perturbation of either the focusing NLS equation, the defocusing NLS equation, or a mixture of both, 
depending on the nature of the initial data. Having already established in §l.l| a connection between MNLS 
and NLS at the level of the corresponding dispersionless systems (JsJ) and (|10[), we want to indicate how this 
connection (and especially the crucial role played by the sign of u(x,t) — c(x,t)) is manifested at the level 
of the inverse-spectral transform theory that underpins both the NLS equation ^ and the MNLS equation 
([I]) as completely integrable partial differential equations. This is a problem of asymptotic analysis in the 
semiclassical limit e -C 1 of the direct (computation of scattering data) and inverse (solution of the matrix 



Riemann-Hilbert problem) components of the inverse-spectral transform algorithm described in 1 1.2 



In Sj2]we present some new asymptotic calculations of the reflection coefficient corresponding to wavepacket 
initial data for the MNLS equation. Then, in §3] and fjljwe interpret these results in the context of, respec- 
tively, the conditions of initial global hyperbolicity and ellipticity of the dispersionless MNLS system (10 1. 



We will show how the type of the latter system at t — influences the asymptotic analysis of the Riemann- 
Hilbert problem of inverse scattering, drawing a comparison with the semiclassical asymptotic analysis of the 
Riemann-Hilbert problems for the defocusing NLS equation and for the focusing NLS equation respectively. 
What will be demonstrated in these sections is quite remarkable: there is a sense in which the direct and 
inverse spectral transforms for the MNLS initial-value problem with globally hyperbolic (respectively ellip- 
tic) initial data are virtually indistinguishable in the semiclassical limit e«l from those for the defocusing 



(respectively focusing) NLS initial- value problem with data obtained from the MNLS initial data via ( 13 1 or 
( 14 ) as appropriate. 

In a final section we will discuss the interesting case of transsonic initial data for the MNLS equation and 
the corresponding phenomenon of a "sonic line" in the semiclassical MNLS dynamics. We also indicate how 
in this situation the asymptotic analysis of the Riemann-Hilbert problem of inverse scattering simultaneously 
contains features characteristic of both the defocusing and focusing cases. 

We are very pleased to be able to contribute our work to this volume of papers in honor of the 85th 
birthday of Peter Lax. Peter's interest in dispersionless limits goes back at least to his consideration of 
von Neumann's conjecture that oscillations produced upon shock formation by a finite difference scheme for 
the integration of a system of partial differential equations of gas dynamics should "on average" (that is, 
in the sense of weak convergence as the grid is refined) reproduce the correct entropy shock profile. In his 
groundbreaking work in the early 1980's with C. D. Levermore on the zero-dispersion limit of the Korteweg- 
de Vries equation [5], Peter showed that von Neumann's conjecture was false (see [5] for more information 
about this interesting problem). 
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2 Semiclassical Analysis of the MNLS Reflection Coefficient 



For initial data <po(x) of wavepacket form Q, the semiclassical parameter e«l enters into the differential 
equation ( 15 ) in two essential ways: as a factor multiplying the derivatives with respect to x on the left- 
hand side, and also in oscillatory exponential factors e ±tS ^ x ^ e residing in the off-diagonal elements of the 
coefficient matrix L on the right-hand side. In this section we present some calculations indicating how the 
limit e — > influences the reflection coefficient s(z;e) corresponding to wavepacket initial data Q appearing 
in the Riemann-Hilbert problem of inverse scattering. 

Let 4>o{x) = (f> t (x,0) have wavepacket form Q. We assume that the corresponding functions po(x) and 
uo(x) defined by ([5| are such that po and u' are Schwartz-class functions, and that po(x) > for all ieE. 
In particular, this condition implies the existence of the limits 



u± 



lim uq(x). 



Letting So '■— 5(0), the relations ^ then imply that 



S(x) = So + / Uo(y) dy = u±x + S± + o(l), i->±oo, 



where 



2.1 Basic setup. 



S+ := So 



±oo 



[uo(y) - u±] dy. 



(36) 



(37) 



(38) 



We will require the Jost matrix solutions for positive real k and positive imaginary k only. Also, it is 
convenient to remove the rapidly oscillating exponential factors e ±lS ( x )/ e from the coefficient matrix L in 
the differential equation (151, so we define related matrices W_|-(j:;z,e) for real z by 

W±(a;; z, e) := exp 



S(x) 
-l— — cr 3 



2c 



J±(x;y/z, e), 
J± (x; i\J—z, e) 



z>0 
z<0, 



(39) 



Thus, the real z-axis corresponds to positive real k (for z > 0) and positive imaginary k (for z < 0). The 
exponential prefactor ensures that the matrices ~W±(x; z, e) satisfy a modification of the differential equation 

"' W± ,MW ±! 



2ae- 



dx 



(40) 



in which the coefficient matrix M = M(a;; z) is independent of e and is defined by 



— 4z + 1 — auo(x) \u^fz\j po{x) 

1 + auo(x) 



Az- 



M := { 



Aa^/zy' po(x) 
— 4z + 1 — auo(x) Ma\J—ZyJ po(x) 



z>0 



z < 0. 



(41) 



Aia\/~z^J po(x) 4z — 1 + auo(x) 
The matrices W±(x; z, e) are then uniquely determined by the normalization conditions 

W-t (x; z, e) exp 



(4z — 1 + au±)x + aS± 

l ; 0"3 



2ae 



±oo, 



z e 



(42) 
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2.2 WKB formalism. 



The Wenzel-Kramers-Brillouin (WKB) method of approximation can be applied to the differential equation 



(40) in the asymptotic limit e J, 0. The form of the WKB ansatz for a column- vector solution w = w(x; z, e) 



ofpO] is 

w = v(x; z, e )e f[x ' z)/(2at \ v(x; z, e) ~ vo(x; 2) + evi(a; z) + e 2 v 2 (x; z) + ■ ■ • (43) 
with the scalar exponent / being independent of e, and the vector v(x] z, e) being expanded in an asymptotic 



scries in the parameter e. Inserting the ansatz into the differential equation (40) and separating powers of e 
results in a hierarchy of equations governing the exponent / and the vector coefficients Vk{x; z): at leading 
order we find 

df 



and then at subsequent orders 



Mv = -i^vo (44) 



Mv fc + i^Vk = -2ia dV f k > 0. (45) 
dx ax 



Clearly, solving (44) with v nowhere vanishing requires that for each real x and z, the characteristic equation 

= (4z - 1 + au (x)) 2 + 16a 2 zp (x) (46) 

ax J 

hold true. We must now distinguish two cases: 

• The oscillatory case. This corresponds to open intervals of x (depending on z £ M) in which the 



right-hand side of (46) is strictly positive. In this case, df/dx will be taken to be one of two purely 



imaginary smooth functions of x: we write df/dx = ±iu>(x; z) where uj(x; z) G R. 

The exponential case. This corresponds to complementary open intervals of x in which the right-hand 



side of (46) is strictly negative. In this case, df/dx will be taken to be one of two purely real smooth 



functions of x: we write df/dx = ±7(2;; z) where 7(2;; z) G R. 

Since po(x) vanishes for large if Az — 1 + an_ 7^ then the oscillatory case holds for x sufficiently 
negative, while if 4z — 1 + au+ ^ then the oscillatory case holds for x sufficiently positive. The isolated 
values of x (depending on z £ M) for which neither the oscillatory nor the exponential case applies because 



the right-hand side of (|46j) vanishes are called turning points. The WKB ansatz ( 43 1 fails near each turning 
point. 

In each oscillatory or exponential interval the derivative df/dx is therefore determined as a smooth 
function of x, and then (44) further requires that Vp reside in the nullspace of the singular matrix M + idf /dx. 



Since the eigenvalues of M are distinct, the matrix M + idf/dx has rank one; let y(x; z) denote a nonzero 
vector in the nullspace chosen arbitrarily, but depending smoothly on x and z. We then see that Vo = voy, 
where vq = vq(x; z) is a scalar factor with smooth dependence on x and z to be determined. Inserting this 



expression on the right-hand side of ( 45 1 for k = 1 shows that dvg /dx will need to be determined as a 
solvability condition for vj.. With dvo/dx so chosen, one then obtains Vi up to the addition of an arbitrary 
nullvector of M+ idf/dx, which may be written in the form viy involving a new scalar factor v±, and dvi/dx 
will then need to be chosen as a solvability condition for v 2 . The procedure continues in this systematic way 
for all k > 0. Clearly, at each order there is produced one new integration constant (independent of x and 
e, but z-dependence is allowed). But these constants may all be combined together, amounting to adding to 
v(a;; z, e) a vector of the form c(z, e)y(x; z) where the scalar c has an asymptotic expansion in powers of e. 

2.3 The reflection coefficient in the case of two real turning points. Barrier 
tunneling problem. Generalization for more than two turning points. 

Let Wjl' and w^' denote the first and second columns, respectively, of the normalized matrix solutions W± 



of (|40|). Since and are linearly independent solutions of (|40|) it follows that there are scattering 
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coefficients T\2{z, e) and T22(z,e) such that the identity 

wf (x; z, e) = T 12 (z, e)w W {x; z, e) + T 22 (z, e)wi 2) (a;; z, e) (47) 
holds. Note that Tjk(z,e) := S'jfc(e" r ' 4 (— iz) 1 / 2 , e) for z 6 R, where the scattering matrix S(fc,e) is defined 



by (18). To compute Ti2(z,e) and T 22 (z, e) we now assume that z € K is such that 

• 4z - 1 + au_ 7^ 0, 

• 4z — 1 + au+ 0, and 

• There exist exactly two turning points X—(z) < and the WKB approximation is in the expo- 
nential case for X—(z) < x < x + (z). Note that this condition implies in particular that z < 0. 

This situation is analogous to that occurring in one-dimensional quantum scattering problems, where a 
particle represented by a wave propagating to the right from x — — oo encounters a potential barrier of 
sufficient height that it lacks the classical energy to penetrate. In the semiclassical limit e <C 1, the wave 
is nearly totally reflected, but an exponentially small fraction of the incoming wave is transmitted through 
quantum tunneling, and appears at x = +oo as a small rightward propagating wave. The mathematical 
problem that arises there is also one of calculating the both the phase of the reflected wave and also the 
exponentially small magnitude of the transmitted wave, and so we refer to the calculation we are about to 



embark upon in the context of the MNLS scattering problem ( 15 ) a barrier tunneling problem. 

(2) 

Our argument is the following: (x;z, e) is approximated by the WKB method (oscillatory case) for x > 

x + (z). If we consider x just to the right of x+(z), we see that w y + '(x; z, e) remains bounded as e 4- 0, although 
the phase varies rapidly. Now, the WKB method fails in a neighborhood of the turning point x = x + (z), 
and to pass through the turning point requires local analysis. However we will be able to proceed with the 
formal asymptotics at leading order simply by assuming that the connection problem is solved by constants 
that are independent of e (as is the case for simple turning points of the equation —e 2 y"(x) + V{x)y{x) = 0, 

(2) 

for example). Under this assumption, the solution w\ (x;z, e) is well-approximated for small e > at x 
just to the left of x+(z) by a linear combination, with coefficients having magnitudes independent of e, of 
the two WKB formulae valid in the exponential region X—(z) < x < x + {z). Continuing this approximation 
to the left, away from the turning point x+(z), we observe that one of the two WKB formulae becomes 
exponentially large compared with the other; therefore if x is any point in the interior of the exponential 
region, wj 2 ' (x\ z, e) is given by 



w + 



/ 1 / +{Z) 

+ (x;z,e) = (y e xp(x;z) + 0(e))exp - — / j(y;z)dy 




(48) 



exponentially small, e i 0, < x < x + (z), 



where 



7(2;; z) := \J — I6a 2 zp (x) — (4z — 1 + au (x)) 2 > 0, %-(z) < x < x + {z) (49) 

and where y CX p(:£; z) is a smooth, e-independent eigenvector of M with eigenvalue i~f(x; z). 

To approximate T 12 (z,e) and T 2 2(z, e) for small e > 0, we will now construct WKB approximations for 
the columns of W_(x;z,e) valid in the oscillatory region x < X-(z). Let 

°~± '■= sgn(4z — 1 + au±), (50) 

and define a real frequency of oscillation by 

u>(x; z) := cr_ \J lQa 2 zpo(x) + (4z — 1 + auo(x)) 2 , x < X-(z). (51) 

For each appropriate value of z € K, this quantity has the fixed sign cr_ throughout its domain of definition, 
and oj(x; z) — > W-(z) as x — > —00, where 

uj-{z) := Az - 1 + «m_ . (52) 
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Finally, let y^ sc (x; z) denote smooth eigenvectors of M(ie; z) with eigenvalues ±lj(x; z) satisfying 



lim yt sc { x -, z ) = 



and lim y osc (x;z) = 



X » X 



(53) 



Then, the leading-order WKB approximations of the columns of W_(x; z, e) valid for x < X— (z) are 

rX-(z) \\ 



w 



(i) 



(x;z,e)=(y- c (x;z) + 0(e))exp(~(cW(z) + ± 



(x;z,e) = (y+ sc (x;z) + 0( £ ))exp * C^\z) 



2a 
1 

2a 



u)(y;z)dy 



x-(z) 



(54) 



uj(y;z)dy 



where (z) and C^ 2 ^ (z) are real constants (independent of x and e) to be determined so that these formulae 
are consistent to order 0(e) with the normalization condition (42 1. We therefore obtain: 



C i - 1 \z) = -C^\z) 



1 

"2a" 



:-(*) i i 

[u}(y;z] -ui-(z)] dy - —uj-(z)x^(z) - -S-. 
„ la I 



(55) 



Now we again invoke the assumption that the constants involved in solving the connection problem across 

(2) 

the turning point x — x_ (z) are independent of e at leading order. Thus, the solution (x; z, e) given just 
to the right of X-(z) by the exponentially large expression (48 1 will also be exponentially large of the same 
magnitude for x just to the left of X-(z), and will be asymptotically represented by a linear combination, 
with constants of the same exponentially large magnitude, of the two oscillatory WKB formulae ( 54 ) valid 
for x < X-{z). This line of reasoning leads to the formulae 



T 12 {z,e) = exp 
T 2 2{z,e) = exp 



2ae 
2a~e 



x+0) 

x-{z) 
x+(z) 



1 (y;z)dy)e^ c{1) ^/t. (l), e\0 



(56) 



-(*) 



j(y;z)dy e"^ 2 '^ • 0(1), eiO. 



The reflection coefficient s(z,e) relevant for the inverse-scattering problem is defined in terms of <Si 2 (fc, e) 
and S22(k,e) by (19) and (22). From the above asymptotic formulae for T 12 and T 2 2 it then follows that 

(57) 



$(z) :=lim(-ielog(s(z,e))) = C (2) (z) - C (1) (z) 



x-(z) ^ 

\oj{y; z) — ijj- (z)} dy H — o;_ (z)x- (z) + S- 



An exact identity within the scattering theory that is a consequence of antiholomorphic spectral symmetry 
and unimodularity of the scattering matrix is 



\T 22 (z,e)\ 2 -\T 12 (z,e)\ 2 = l, z < 0, 
and this together with the asymptotic formulae for Ti 2 and T 22 implies that 



t(z) :=-hm(elog(l-| S (z, e )| 2 )) 



x + (z) 



x-{z) 



i(y; z) dy 



(58) 



(59) 



whenever z is a (necessarily negative) value of the spectral parameter for which there are exactly two real 
turning points. To summarize, for such z < 0, the reflection coefficient s(z, e) has modulus exponentially 
close to 1 (as measured by r(z) > 0) and real fast phase approximated by e~ 1 ^(z). 
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The functions &(z) and t(z) can be related by analytic continuation in the case that the coefficient 
functions uq and po are analytic. Indeed, suppose this is the case, and that zq is a real value of z at which 
the real turning points x±(z) coalesce from below (z < zq), so that in particular x+(zq) = X-(zo). For 
z < zq, x±(z) are two branches of the inverse of an analytic function z = z{x) satisfying z' (x±(zq)) = 0. If 
we presume the generic condition z" (x±(zq)) < 0, then it follows that x±(z) behave locally like square roots 
of zq — z. In particular, a full circuit of z about zq in the positive (negative) sense results in a half-circuit 
of x±(z) about x±(zo) in the same sense, and this also leads to the permutation of the two turning points: 
x±(zq + (z — zo)e 27Tt ) = x T (z) holds for z < zo, as does x±(zq + (z — zo)e~ 27 ™) = x T (z). For z < zq, we 
choose xq < X- (z) and write $(2) in the form 

$(z) = - / [u(y;z) - co^(z)]dy+ -x uj_(z) + + - / u(y;z)dy, z < z . (60) 
a J -00 " " Jx 

The first two terms on the right-hand side are analytic in z at z = zq, and hence will return to the same 
value upon a complete circuit of z about Zq in either the positive or negative sense. The third term, however, 
experiences monodromy upon such an analytic continuation, because the turning points are exchanged: 

^ rx„(z + (z-zo)e ±27rl ) ^ /-x_(z) ^ r x +i z ) 

u{y; z + (z - z )e ±2vl ) dy = ui(y;z)dy+~ uj T (y;z)dy, (61) 



O: 



where for X—(z) < y < x+(z), 



u±(v,z) ■= limu)(y±iS;z) = T^-l(y,z)- (62) 
0^0 



It therefore follows that 

<P{z + (z - z )e ±27rl ) = $(z)±ia_T(z), z<z . (63) 

It is easy to check that exactly the same formula holds true for analytic continuation around a point zq G K 
at which x±(z) coalesce like square roots from above [z > zq). Even simpler reasoning (write r as a loop 
integral around a complete branch cut between X—(z) and x+(z)) produces the continuation formula 

t(z + (z- z )e ±2 ™) = t(z) (64) 

for analytic continuation about both types of branching points zq . 

Finally, we note that the arguments presented in this section also carry over to the case when there exist 
arbitrarily many simple turning points (necessarily an even number). Indeed, if x_(z) is simply reinterpreted 



as the left-most turning point, then the formula (57 1 for the phase $(2) of the reflection coefficient still holds, 



and the formula (59 1 for t(z) needs only to be modified by replacing the integration interval by the union 



of intervals in which the exponential case holds, that is, where the radicand of 7(0;; z) as defined by (49 1 is 
positive. 



2.4 



The reflection coefficient in the absence of real turning points, 
barrier reflection. 



Above- 



z € 



is a value for which the right-hand side of ( 46 1 is uniformly bounded away from 

^ as 



Now suppose that 

zero for x £ M. In fact the right-hand side of ( |46[ ) is necessarily positive in this case, since po(x 
\x 



00. Therefore we are in the oscillatory case for the WKB method over the whole real a;-axis, and in 



particular, for smooth po and uq the function u)(x; z) is well-defined by (51 1 as a smooth function of x G K 
whenever zgl. Moreover, it can be shown that the WKB ansatz is uniformly valid for x G K to all orders 
of accuracy. In other words, for each choice of sign in the formula df/dx = ±iu>(x;z) there exists a one- 



dimensional subspace of true solutions w(x;z,e) for which the asymptotic expansion (431 holds uniformly 
for the magnitude of the difference between v[x\ z, e) := w(a;; z, e)e~^ x ' ,z,t ^ ty2a '^ and the partial sum 

Vg(a;; z) + evi(x; z) + ■ ■ ■ + e N \rpj(x; z) has a maximum value over x G K that is 0(e N+1 ) as e J, 0. 
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The absence of turning points in a quantum scattering problem occurs when a classical particle incident 
on a potential barrier from x = —oo has sufficient energy to penetrate the barrier (or propagate "above" 
the barrier) and arrive at x = +00. Quantum mechanically, however, when the particle is represented as 
a right-going incident wave, the barrier still causes the generation of an exponentially small reflected wave 
(and similarly, the transmitted wave is attenuated by an exponentially small fraction). In the context of the 
scattering problem (15 1 for the MNLS equation, the "reflection coefficient" s(z, e) can correspond either to 
the magnitude of the physical reflection coefficient or to its reciprocal, depending on a monodromy index 
associated with the coefficients representing the barrier. Despite this difference, we refer to the calculation 
we are about to embark upon as a problem of above-barrier reflection. 

Let N± = N± (x; z) be given by 



{y/ '16a 2 z po(x) + [Az — 1 + auo(x) ± ui{x; z)] 2 , z > 

y/± (16a 2 zp (x) + [4z - 1 + au (x) ± uj(x; z)] 2 ), z < 0. 



(65) 



In both cases the square roots are taken to be positive. (It is obvious that the radicand appearing in the 
definition of N± is positive if z > 0. On the other hand, if z < 0, we have 



16a 2 zpo(x) + [4z — 1 + auo{x) ± lj(x; z)] 2 = 2ui(x- 1 z) [ui(x; z) ± (Az — 1 + auo(x))] , 



(66) 



and since z < 0, < \cj(x, z)\ < \4z — 1 + auo(x)\ so the sign of (66) is the same as that of uj(x; z), namely 
er_, times that of ±(4z — 1 + auo(x)). But since we are in the oscillatory case for all x £ K by assumption, 
z < implies that Az — 1 + auo(x) is bounded away from zero for and therefore the latter sign can be 

taken from the limiting value at x = —00, namely sgn(±(4z — 1 + cxuq(x))) = ±cr_, so the radicand is again 
positive.) An eigenvector of M with eigenvalue ±u(x; z) assumed real is given by 



y±(x;z) 



AayfzyJ Po(x) 

Az — 1 + auo(x) ± ui(x; z) 



Az 



Aict\J-~Zy/ po{x) 
- 1 + auo(x) ± ui(x; z) 



z > 



z < 0. 



(67) 



It is then a straightforward calculation to confirm that by choice of the normalization constants N± , we have 



z > 
z < 0. 



y±(x; z) T y T (a;; z) = and y±(z; z) T y ± (x; z) = 
We also record here the asymptotic values of y±(x; z) as \x\ — > 00: if z > 0, we have 



lim y + (x]z) = 



and 



lim y_ (x\ z) — 



while 



x— >-+oo 





0~ 


< 


1 




1 








and 



lim y_ (x; z) = < 

x— ^+00 



(7+ ^ <T_ 



(68) 



(69) 



(70) 



On the other hand, if z < 0, then we necessarily have a + = er_ under the assumption that there are no real 
turning points, and so 

"0" 



lim y+{x;z) 

\x\— »oo 



and 



lim y-{x;z) 

\x\—>oo 



(71) 
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Now let Y(x;z) := (y + (x; z), y_ (x; z)) be the eigenvector matrix. From (68) it follows that 

_! _ |y(x;z) t , z > 



Y{x;zY 



\a 3 Y(x; Z y 



z < 0. 



From the identity Y(x; z) 1 Y(x; z) = I we get by differentiation 

= ' !YiJ r' Vr: z) + Y(x; z)" 1 dY( *' z) 



fix 



so combining with ( 72 ) we see that regardless of whether z > or z < 0, 

_ x aY(x; z) 



Y(a:;z)- 



fix 



h(x; z) 
—h(x;z) 



(72) 



(73) 



(74) 



for some smooth function h(x;z). We would like to use this information to calculate the Jost solution 

(2) \ L 

w+ (x; z, e) in the limit e I 0. According to (42 ), we should choose 

f ix; z) — iuS+ + it7 + tj-U) + (z)x — ia + <7- I [uj(y; z) — cu + (z)] dy to obtain w^(x; z, e). (75) 



Here, 



w+(z) :— lim u;(x; z) = (7 + <7_(4z — 1 + ait+). 

a;— f +00 



(76) 



It is now easy to carry out the calculation of the WKB approximation to all orders, by asymptotically 
diagonalizing the system governing v(x;z,e) := w(x; z, e^ - -^'^/' 2 " 6 ) with the use of the substitution 
v(x; z, e) = Y(x; z)z(x; z, e). In the case <r+ = er_ so that df /dx — iu(x; z), the system satisfied exactly by 
z = (Ci, C2) T consists of the equations 



= — 2iae ( h(x; z)(2(x; z, e) H — r~(x; z, e) 
\ dx 



2~-(.r: : (.)•::.( ) -2/n< j — ft,(x; z)d (x; z, e) + (x; z, e) 



(77) 



On the other hand, in the case cr + 7^ cr_ (implying z > 0) so that df /dx = —iuj{x; z), the equations become 



2w(x; z)£i(x; z, e) = — 2iae ( h(x; z)(2(x; z, e) H — r~fe z, e) 

\ ax 



dC,2 

= — 2iae ( —h(x; z)£i(x; z, e) + -j— (x; z, e) ) . 



(78) 



Note that since p and u' are Schwartz class, for all z under consideration, u>(x; z) _1 is bounded and infinitely 
differentiable, and h(x; z) is Schwartz class. 

To solve (77) we eliminate Ci using the first equation and building in the correct boundary condition 

(79) 



implied by (42) at x = +00: 



Ci(x;z,e) = 1+ / h(y;z)C 2 (y,z,e)dy, 



and therefore rewrite the second equation of ( 77 ) in the form 

~d( 2 



C2(z;z,e) 



icee 
oj(x; z) 



dx 



(x;z,e)-h(x;z)-h(x;z) / h(y; z)( 2 (y; z,e) dy 



-f- OG 



(80) 
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Iterating this equation starting with the initial guess £2 (x; Z, e) = produces the partial sums of the asymp- 
totic power series for £2(2; z; e): 

00 

C 2 (x;z, £ )^^C2,n(x;z)e", e | 0. (81) 

n=l 

It is straightforward to confirm by induction that every term in this power series is a Schwartz-class function 
of x. 



In a completely analogous fashion we solve ( 78 ) obtaining the asymptotic power series 



Ci(x;z,e)^J2c 1}n (x;z)e n , eiO (82) 

n=l 

in which each term is a Schwartz-class function of x. 

The uniform validity of the WKB approximation in the absence of turning points means that we may 
interchange limits to deduce the behavior of w\ (x;z,e) as x — > —00. Now, the coefficients Tyi{z,e) and 
^22 (2, e) in (47) may be expressed in terms of Wronskians as 

T\i{z, e) = det {v*^ (x; z, e), wf 2 -* (x; z, e)^j and X^z, e ) = det ^w^^x; z, e), w+\x; z, e)J (83) 

and the right-hand sides are actually independent of x, so we may take the limit x — > —00 and use the 



asymptotics given by ( 42 ) to write these formulae exactly in the form 

T 12 (z,e)= lim e^-^+^-V^w^lix-z^e) 



X—t — OO 



(84) 

T 22 (z,e) = lim e-^-^ +aS -^ (2a ^w^ 2 (x;z,e). 



x^ — oo 



("2) 

Substituting the WKB asymptotics for the components of w\ 1 (x; z, e) we find differing behavior of s(z, e) 
depending on the sign of the monodromy index m := (T + (T_: 

• If m = +1, (that is, cr + = cr_), then T±2(z, e) = o(e n ) for every integer n while T22(z, e) is bounded 
away from zero in the limit e 1 0. The reflection coefficient s(z,e) = —Ti2(z,e)/T 2 2(z,e) is therefore 
small beyond all orders in e. 

• If m = — 1 or <j_|_ 7^ <r_ (which implies in particular z > in the absence of real turning points), then 
T 2 2{z,e) = o(e") for every integer n while Tx 2 {z, e) is bounded away from zero as e 4- 0. In this case, 
the reciprocal of the reflection coefficient is small beyond all orders in e, and we should expect the 
reflection coefficient to be exponentially large. 

In the case that s(z,e) is small (cr + = cr_), the jump matrix in the Riemann-Hilbert problem of inverse 
scattering is well-approximated by the identity, and no further information about the reflection coefficient 
is really required. However when s(z, e) is large (cr+ ^ <r_), further asymptotic analysis of the Riemann- 
Hilbert problem requires a leading-order formula for the reflection coefficient. This is a challenging problem 
of exponential asymptotics that appears to require the analyticity of po and uq with respect to x as the most 
effective methods rely on analytic continuation of the WKB solutions from the real line into the complex 
x-plane. If po an d uq are analytic functions of x, then given z€K for which there are no real turning points 
we may instead seek complex turning points, that is, complex roots x = x{z) £ C of the right-hand side of 



(46). These obviously come in complex-conjugate pairs for real z. If the right-hand side of (46) is analytic 



in a strip of the complex x-plane containing the pair of complex-conjugate turning points closest to the 



real x-axis, then one can show that the formula (57) continues to hold for the phase $(z) of the reflection 
coefficient, in which x_(z) is interpreted as the turning point in the upper (respectively, lower) half-plane 
for o — = — 1 (respectively, er_ = +1), and the path of integration is confined to the strip of analyticity. It is 
then clear that 3?{i$(z)} > 0, making the reflection coefficient exponentially large. 
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For z real and of sufficiently large absolute value, it is obvious that there exist no real turning points, 
and also that cr+ = er_, making the reflection coefficient small beyond all orders. Also, it is easy to see that 
cr_|_ 7^ <r_ exactly for those z £ K lying in the interval / given by 

I .= {— (1 — a max{u + , «-}), — (1 — a min{w + , w_})^ . (85) 

Note, however, that z € I does not guarantee the absence of real turning points; for example for each z £ I 
with z < there necessarily exists at least one real turning point (in fact at least two if they are simple) . 

2.5 Discrete spectrum. 



Locating eigenvalues for the scattering problem ( 15 ) is also a difficult problem of exponential asymptotics. 



The discrete spectrum (at least the non-embedded part) is nonreal, and it is easy to see from the characteristic 



equation ( 46 ) that in order to have exponential decay of the eigenfunction as x — > ±oo it is necessary to 
have 7^ 0. To construct an eigenfunction using the WKB method, one needs to be in the exponential 

case for sufficiently large \x\ and to connect the two regions of decay through a central oscillatory region via 
two (or possibly more) turning points. The problem is that for typical nonreal z, there exist no real turning 
points at all. One way to get around this difficulty is to further assume that the initial data (po(x), uq(x)) 
are real analytic functions of x admitting analytic continuation into the complex x-planc. In this situation, 
it is frequently possible to find pairs of complex turning points x±(z) £ C and to construct an approximate 
eigenfunction that is oscillatory along a certain curve connecting the turning points and that is exponentially 
decaying away from the turning points in basins that contain the two points at infinity on the real x-axis. 
The condition that a curve exists in the complex x-plane connecting the turning points along which the 
real part of the WKB exponent is not changing determines a curve in the complex z-plane near which one 
expects to find 0(e _1 ) eigenvalues approximately determined by a Bohr-Sommerfeld integral condition that 
quantizes the phase increment between the turning points. This sort of calculation is described in the context 
of the nonselfadjoint Zakharov-Shabat eigenvalue problem in |11) . 

If the reflection coefficient s(z,e) = —Ti2(z,e)/T22(z,e) has a meromorphic continuation into the upper 
half z-plane from the neighborhood of some zo £ K (which is really a statement about the scattering coefficient 
T\2(z, e), since 122(2, e) is analytic in the upper half-plane), then it will have simple poles at the eigenvalues 



= kj, and the corresponding proportionality constant jj defined by (201 is given by 7j = Ti 2 (zj,e). This 



observation can allow the direct computation of eigenvalues to be avoided in some situations. 



3 Supersonic Initial Data 

Here we consider the solution of the MNLS equation by the inverse-spectral transform in the case that 
the initial data is globally supersonic for the dispersionless limit system (10 1. That is, we assume the initial 
data to satisfy the condition that Q(x) := a 2 po(x) + auo(x) — 1 > for all x £ M. 

The discriminant of the quadratic polynomial in z on the right-hand side of the characteristic equation 



(46) is 256a 2 po(x)Q(x). Hence for supersonic initial data there exist two real roots z = z±(x) defined for 
all x £ K, with z + {x) > Z—(x). Also, directly from the quadratic formula one checks that the supersonic 
condition and the positivity of po(x) imply that z+(x) < for all x £ K, with equality only for those x 
for which auo(x) = 1. In general, the locus of complex values of z parametrized by x £ K for which the 



right-hand side of the characteristic equation ( 46 1 vanishes is called the turning point curve. So we see that 
for globally supersonic initial data the turning point curve is more specifically a real curve. 

Due to the presence of the additional factor of po(x) in the discriminant, we also have that z + {x) and 
Z-(x) coalesce as x — > ±00 to the value |(1 — au±). Therefore, upon setting 

zj := ramz-(x) and z-o := maxz+(x) < 0, (86) 

we have the following: 
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• For zl < z < zr , there exist real turning points (and by Sard's Theorem the turning points are simple 
for z in a subset of this interval of full measure). For such z the reflection coefficient s(z, e) has modulus 
exponentially close to 1 as measured by the function t(z) > defined by (59 1, and has approximate 
real-valued phase e _1 $(z) where <&(z) is defined by (57). 



For z < zl and z > zr, there exist no real turning points at all, and since zl < t(1 — au±) < zr, we 



have / C [zl, zr] where / is the interval defined by (85 ). Consequently, the reflection coefficient s(z, e) 
is small beyond all orders for z [zt, , zr] . 

Moreover, one of the general results of [2] is a hard estimate on the discrete spectrum of the scattering 



problem (15) that confines the eigenvalues to the "shadow" cast by the turning point curve in the complex 



z-plane by light projected from infinity along vertical lines toward the real axis, up to an error proportional 
to e. In the present case, the turning point curve itself is real, and hence it is its own shadow. As this 
shadow has no intersection with the open upper half-plane, the conclusion is that, at least at the level of 
the semiclassical approximation, there are no eigenvalues at all for supersonic initial conditions of the type 
considered here. 

These results should be compared with those for the direct scattering problem for the inverse-spectral 
transform adapted to the defocusing NLS equation (pi) with k = —1). The analogue of (15) in this case is 



the self-adjoint eigenvalue problem (27) in the case n = — 1. The continuous spectrum for this problem is 



the real axis of the spectral parameter (in this case, A), just as in the case of the scattering problem ( 15 1 for 
MNLS. By self-adjointness, there can be no nonreal eigenvalues (discrete spectrum) whatsoever, which is an 



exact version of the semiclassical asymptotic "shadow" estimate for the MNLS scattering problem (15) with 
supersonic initial data. 

Our expectation is that the semiclassical defocusing NLS dynamics, after being written in terms of the 



MNLS fields p e and u e via the explicit mapping (13), should resemble the semiclassical MNLS dynamics in 



the case of supersonic initial data. In other words, we begin with functions po(x) and Uq(x) satisfying the 
condition Q(x) > for all x £ R, and consider at the same time the solution (f) t (x,f) of the initial- value 
problem for the MNLS equation ([I]) with wavepacket initial data given by Q and and the solution 
4> t (x,t) of the initial-value problem for the defocusing NLS equation ([2]) (with k = —1) with wavepacket 
initial data 4>o(%) given by 



Mx) = A(x)e i§ ^, A(x) := y/po(x)Q(x), S(x) := S + / [u Q (y) + 2ap (y)} dy 



(87) 



For convenience, we choose the integration constant So so that S(x) — S(x) + o(l) as x — > — oo, that is, we 
set 



So := So + 2a 



Po(y)dy. 



We have already described the scattering data for the MNLS direct spectral problem ( 15 ) with potential 



0o ; and now we must do the same for the Zakharov-Shabat spectral problem (27) in the self-adjoint case 
of k = — 1, with the potential function 4>q. We begin by introducing a shift and a-dependent scaling in the 
Zakharov-Shabat spectral parameter by writing 



X = X(z):= 2 -(z- 1 - 
To prepare this spectral problem for WKB analysis we make the transformation 



(89) 



W 



zs,-i 
± 



(x, z, 



exp 



.S(x) 



j 



zs,- 



1 (x;X(z),e) 



(90) 



analogous to (39) in (27) for k = — 1 (subject to the above substitutions), resulting in the differential equation 



2ae- 



dW 



zs. 
± 

dx 



= iM^-^wf' 



(91) 
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with coefficient matrix 



M 



zs,-i 



— 4z + 1 — au Q (x) — 2a 2 p (x) 
-2ia^f pu(x)Q(x) 



-2iay/p (x)Q(x) 
4z — 1 + auo(x) + 2a 2 pq(x) 



(92) 



Making the WKB ansatz for coiumn vector solutions of the form w = v(x; z, e ) e f(x;z)/(^) l ea ds in the same 
way as in <|2]to exactly the same characteristic equation (46) as governs f in the MNLS case. In other words, 
M zs ' _1 defined by (92) and M defined by (41) are similar matrices (in fact they become equal in the limits 
x — > ±oo). 

By following the line of reasoning described in detail for the MNLS spectral problem in S|2j one sees easily 
that the quantity 



,zs,- 



„ZS,-1 



(X(z),e), 



z e 



(93) 



is small beyond all orders in e for z \z\,, zr], while for z € (zl, 2r) we have the approximations s ZS ' x (z, e) 



e »*(*)/e and i 



e zs : - 



-rW/e 



where <&(z) and t(z) are given by exactly the same formulae 



((57) and (59)) as in the MNLS case and where [#l,Zr] is the interval in which there exist real turning 
points. In other words, at the level of semiclassical approximation, the MNLS reflection coefficient s(z; e) is 
indistinguishable from the defocusing NLS reflection coefficient s ' _1 (z;e). 

For one thing, these considerations show that for supersonic initial data {po,Uo) used to construct 
wavepacket initial data for the MNLS equation, the Riemann-Hilbert problem of inverse scattering as de- 



scribed in 1 1.2.1 simplifies considerably. Neglecting the discrete spectrum means that there are no poles in 
the matrix unknown N(z). Moreover, for z < outside of the interval [zl, Zr] the jump condition on the real 
axis simplifies to N+(z) = i cr3 N_(z)« _<T3 , while for z > outside of the interval [zl, zrJ the jump condition 
on the real axis simplifies to the trivial jump N + (z) = N_(z). Only in the negative interval [zl,zr] is the 
jump condition nontrivial, taking the approximate form 



N+(z) = i CT3 N_(z)r 



I _ e i(20(z:x,t) + *( Z ))/f 

e -»(2fl( Z ;x,i)+*(»))/e e -r(z)/e 



z e (z L , zr). 



(94) 



However, equally impressive is the analogy with the Riemann-Hilbert problem for defocusing NLS with 
wavepacket initial data of the form (87). Indeed, due to the identity zs (A(z); x, t) = 8(z;x,t), the jump 



condition for the Riemann-Hilbert problem for the latter equation as described in ^1.2.2 (and written in 



terms of the complex variable z rather than A) is negligible for z ^ [zl, zr] while for z G (zl, Zr) one has the 
approximate jump 



N+ S ' _1 (A(z)) = N zs ^ 1 (A(z)) 



1 

-i(29(z;x,t)+&(z))/e 



J(2e{z;x,t)+$(z))/e 
-r(z)/e 



z e (z L ,z R ). 



(95) 



In both jump conditions (94) and (95) there are correction terms we have not written on the off-diagonal 



that come into play near the endpoints zl and zr to bring the jump matrix smoothly to the identity; note 
that t(zl) = t(zr) = due to coalescence of turning points. This is a technicality to be handled by the 
installation of appropriate parametrices in the complex plane near these two points and is not expected 
to affect the results at leading order in the semiclassical limit. Setting this technicality aside, we notice 
that the "limiting" Riemann-Hilbert problems for MNLS and defocusing NLS are nearly identical, with the 
differences being: 

• The point of normalization to the identity for the MNLS Riemann-Hilbert problem is z = 0, while that 
for the defocusing NLS Riemann-Hilbert problem is z = oo. 

• The jump condition for the MNLS Riemann-Hilbert problem involves the additional effect of conjuga- 
tion of N_ (z) by i a3 for all z < 0. 

Regardless of these differences, the same basic asymptotic technique can be applied to both of these problems. 
The key technique first appeared in the study of Deift, Venakides, and Zhou [T] of the small e behavior of 
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the solution of the Cauchy problem for the Korteweg-de Vries equation 



du e 
~dt 



+ 



dx 



1 , 

-K 

2 e 



dx 3 



u £ (x,0) = u (x), 



(96) 



another kind of problem whose e = limit equation is of hyperbolic type, and for which the associated linear 
eigenvalue problem in the inverse-spectral transform is self-adjoint. We first make parallel substitutions of 
the form 0(z) := N(z)e i9 ^ a3 / e an d O zs ~ l {z) := N zs - 1 (\{z))e l9 ^ a3 ' t for a suitable and common scalar 
function g{z) analytic for z G C \ [zl, Zr] and bounded at z = oo. Since the normalization conditions differ 
for the MNLS and defocusing NLS problems, we choose (arbitrarily) to preserve the normalization condition 
in the MNLS case by imposing the condition that g(0) = 0. It is important for semiclassical analysis 
that g is to be chosen independent of e. Now, this substitution preserves the domain of analyticity of the 
unknown matrix functions, and also the MNLS normalization: 0(0) = I. The defocusing NLS normalization 
condition becomes O zs,_1 (cx)) = e 19 ^ 00 ^ 3 ^ . The substitution also preserves the trivial jump conditions 
0+(z) = 0-(z) for z > 0, 0+(z) = i CT3 0_(z)iT CT3 for z < with z <£ [z L ,z R ], and 0+ S = ZS 1-1 (z) 
for z G" [zl, zr]. The jump conditions for z G (zl, zr) take a similar form: 



0+(z) = i CTs O_(z)r 



and 



where 



O 



zs ' _1 (z) = O^'-^z) 



3 »(g+(*)-fl_(«))/e 
-i(Ji + (*)+M*))A 



e *(9+( 2; )-5-( 2 ))/ e 
,-*()i + (*)+M*))A 



_ e i(ft + ( Z )+h_(z))/e 
,,i(iT(»)-(ff + (*)-g_(a)))/« 



_ e »(fc + (*)+M*))A 

„t(iT(*)-(fl + (z)-fl_ («)))/« 



z € (z L , z R ), 



z e (z L ,z R ), 



h(z) := 6{z-x,t) + -<f>(z) - g(z), z G C\ [z L ,z R ]. 



(97) 



(98) 



(99) 



We insist that g(z*) = g(z)* , which makes g+{z) — g~(z) imaginary and h+(z) + h—(z) real. Now the idea 
described in pQ is to try to choose g{z) so that (zl,zr) is partitioned into a finite union of three types of 
subintervals: 



• "Voids," subintervals in which g^ 

• "Bands ," subintervals in which h 



g_ = and 1%+ + h— = 2h is strictly increasing. 
_l_ -r h- is constant and < —i(g+ — .<?-) < r. 



• "Saturated regions," subintervals in which g + — g_ = it and h + + h- is strictly decreasing. 

Via appropriate factorizations of the jump matrix in voids and saturated regions and subsequent "steepest 
descent" deformations in thin "lenses" about these two types of intervals pQ it can be seen that these intervals 
do not contribute at leading order to the solution and can be neglected (again, modulo parametrices of Airy 



type at band endpoints). In each band subinterval however, the diagonal elements of the jump matrix in (97) 
and (98) decay exponentially to zero as e — > 0, while h + + h- is equal to a fixed real constant (independent 
of z, but generally depending on x and t). Let h + + h- = ( n G M for z e (<z„, b n ) C (zl,Zr), the n th of N 
band subintervals, counted from left to right. This reasoning leads to two related model Riemann-Hilbert 
problems for the MNLS and defocusing NLS cases. The model problem for the MNLS case is to find a 2 x 2 
matrix O(z) with the following properties: 

Analyticity: O(z) is analytic for z € C \ (— oo,0] and takes boundary values that are continuous 
except at the band endpoints {amb n }n=o wnere the matrix elements of O(z) are allowed to blow up 
like at worst a negative one-fourth power singularity. 



Jump condition: For z < outside all bands we have 0+(z) = i°" 3 0_(z)i <T3 . 
instead: 



In the bands we have 



0+(z) 



5 6_(z)r CT3 



o 

-iCn/e 



a n < z < b T , 



o,. 



,N-1. 



(100) 
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Normalization: 0(0) = I. 

On the other hand, the model problem for the defocusing NLS case is to find O zs ^ 1 (z) characterized by 
the following conditions: 

Analyticity: O zs,_1 (z) is analytic for z € C \ U^T 1 [a n ,6 n ] and takes boundary values that are 



continuous except at the band cndpoints {a ni b n } n=0 where the matrix elements of 
allowed to blow up like at worst a negative one-fourth power singularity. 

Jump condition: 

_ e iC»/ 



»zs,-i, 



6f'- 1 {z) = 6 zs >- 1 (z) 



o 

e -iCn/e 







a n < z < b n , n = 0, . . . , N — 1. 



z arc 



(101) 



Normalization: O 28 '- 1 ^) = e »9(°°)^A. 



These two model problems are explicitly and uniquely solvable in terms of Riemann theta functions of genus 
N — 1. Let us describe this procedure in detail in the simplest case, N = 1. In this case, the solution of the 
defocusing NLS model problem is completely standard: 



O 28 '" 1 ^) 



e *ff(oo)CT 3 /e e iCoo-3/(2e) 



1 1 

i —i 



(z - boY*l\z - a )" CT3 



/4 



1 1 

z — i 



-iCo<T3/(2e) 



(102) 



where the principal branch of the one-fourth power functions is taken in each case. The solution of the MNLS 
model problem is somewhat nonstandard due to the conjugation by i 173 of the boundary value 6 zs,_1 (z). It 
is perhaps easiest to solve this one by reverting from the z-plane to the fc-plane, where k = z 1 ! 2 (principal 
branch). Indeed, if we define 

'6(fc 2 ), 3?{fc}>0 
i a3 6(k 2 )i-' 73 1 3£{fc}<0, 



6(k) := 



(103) 



then d(k) has two cuts, one on the positive imaginary fc-axis between the points i\f—bo and i^—ao, and the 
other on the negative imaginary fc-axis between — is/— ao and —i\/—bo (note that ao < bo < 0). If we denote 
the boundary value taken on the cuts from the left (respectively right) half-plane as + (fc) (respectively 
0_(fc)), then in both cuts the jump condition turns out to be exactly the same: 



6+(fc) = 6_(fc) 





_ e -iCo/e 







(104) 



and of course we have the normalization condition 0(0) = I. Solving for O(fc) with jumps of this form is 
now a standard procedure, and the solution is 



6(A) 



1 1 

i —i 



H(fc - i^b )y*'\-i(k + i^a7o)Y 3/i 

■ (~i(k + *v/^Mr ff3/4 H(fc - iV^ao)r a3/4 



i i 

i —i 



where again the principal branch is intended for all power functions (in particular, this makes the cuts lie 
on the imaginary axis as intended). To obtain O(z) we simply restrict k to the right half-plane and write 
k = z 1 / 2 (principal branch). Thus, 



O(z) 



= iCo<T 3 /(2e) 



1 1 

i —i 



(-liz 1 ' 2 - iv /=^))-3/4 ( _^l/2 + iy /=^ys/4 
{-liz 1 ' 2 + z v ^))-^/ 4 (-z(z 1 / 2 - z^))-^ 4 



1 1 

i —i 



-iC a 3 /(2e) 



. (106) 
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A completely rigorous semiclassical analysis comes from comparing global parametrices for O(z) and O zs, ^ 1 (2:) 
with the unknown matrices themselves and proving that the discrepancy in each case is the solution of a 
"small norm" Riemann-Hilbert problem, that is, one with near-identity jump matrices. In each case, the 
global parametrix is exactly equal for \z\ sufficiently large to the solution of the corresponding model prob- 
lem we have just presented for N = 1. This fact allows us to extract approximate formulae for cf> e (x,t) and 
<j) t (x,t) that are valid in the semiclassical limit e -C 1, by means of the formulae (26) and (35). It only 
remains to multiply the model solution on the right by e _I9 ^' CT3 / e and calculate asymptotics as z — > oo. In 
the defocusing NLS case we have 



4> e (x,t) ~ 2i lim A(z) (6 zs ^ 1 (z)e' t9{z),T3/e ) 

2— S-OO \ / ] 



bo - a 



,i(2 9 (oo)+Co)A 



12 a 

from which we obtain approximations to the corresponding Madelung variables in the form: 

9 



p e (x, t) ~ p(x, t) := — - — „ a ° and u € (x,t) ~ u{x,t) 



On the other hand, in the MNLS case we have 



dx 



(2 5 (oo) + Co) 



- lim z 1/2 



6{z)e 



-ig(z)a 3 /t 



(107) 



(108) 



(109) 



from which we obtain the Madelung variables 



Pe(x,t) ~ p(x,t) 



and u e (x, t) ~ u(x, t) :- 



dx ' 



(110) 



Now, the endpoints (a , bo) and the phase constant £ are certain functions of x and t as required to construct 
the function g. Indeed, in the case N = 1 the z-derivative of g is necessarily given by the following formula: 



S(z) 
2ni 



b 



<f>'(s)ds 
S + (s){s - z) 



it'(s) ds 

saturated regions S(s)(s — z) 



+ 9'(z;x,t) + 



8tS(z) 



(111) 



where S(z) 2 = [z — ao)[z — bo), S(z) has its branch cut in the interval [ao,6o], and S(z) = z + 0(1) as 
z — > oo. The second integral could be not present at all, or could involve one or both of the intervals 
(zL,a ) and (6 ,zr). For arbitrary (ao,&o)> this formula has a Laurent expansion for large z of the form 
g'(z) — M (a ,bo;x,t) + M\{aa, bo] x, t)z~ l + 0(z~ 2 ). The condition that g'{z) = 0(z~ 2 ) as z — > oo 
(necessary to ensure that g(oo) is finite) then implies that (ao, &o) should be determined as functions of (x,t) 
via the moment conditions 



Mo(a,Q, bo; x, t) = and M^ag, bo; x, t) = 0. 



(112) 



With the correct configuration of voids and saturated regions, the band endpoints are determined from these 
conditions as smooth functions of (x, t) by continuation (via the implicit function theorem) from t = 0, at 
which time S(z) 2 is taken to be the polynomial 



S(z) 2 = — \Uz - 1 + auo(x)) 2 + I6a 2 zp (x)] , t = 0. 
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(113) 



(In other words, at t — 0, the points ao(x, 0) and bo(x, 0) trace out the two branches of the real turning point 
curve.) Now with the endpoints determined as functions of (x,t), it is easy to write down a formula for the 
x-derivative of g(z) (differentiation with respect to x removes the complicated functions $(z) and r(z) from 
the jump conditions): 

: g(z) = —(z;x,t)--^ + -S(z). (114) 



9a;'" 



dx 



2 dx 
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Since g(0) = independently of x, we must impose 



dx 



(3(0)) = which implies 



(115) 



and then it follows easily that 



d_ 
dx 



5(00) 



(116) 



The equations (115) and (116 1 then prove that the semiclassical approximations for p t (x,t) and u £ (x,t) 
defined by (108) and those for p e (x,t) and u e (x, t) defined by (110) are linked by the identities (13) for 
all (x, t) for which the single-band ^-function serves to asymptotically and simultaneously reduce both 
Riemann-Hilbcrt problems (a condition requiring the auxiliary inequalities on the boundary values of g to 
be confirmed). Moreover, the four functions p(x,t), u(x,t) 7 p(x,t), and u(x,t) are obviously explicitly given 
in terms of ao(x, t) and bo(x, t), the latter being determined by a smooth solution to the (algebraic) moment 



conditions connecting to given supersonic MNLS initial data by the identity (113). It is an exercise to 



confirm by differentiation of the moment conditions with respect to x and t that p(x,t) and u(x,t) satisfy 
the dispersionless MNLS system (10), while p(x,t) and u(x,t) satisfy the dispersionless NLS system ^ in 
the defocusing case of k = — 1. In other words, the moment conditions provide the solution of the Cauchy 
initial- value problem for both (l8|) and ( 10 1 in implicit — but algebraic — form. 



The significance of the analogy between the Riemann-Hilbert problem for MNLS with supersonic initial 
data and that for defocusing NLS with corresponding wavepacket initial data goes far beyond the case N = 1 
however. Indeed, the fact that the very same g-function serves to simultaneously reduce both Riemann-Hilbert 
problems in the semiclassical limit implies immediately that, for example, when the (essentially common, 
via the mapping (13)) solution to the hyperbolic nonlinear dispersionless systems ^ (for k = —1) and 



(10) forms shocks that are then dispersively regularized by the neglected term d x F[p], the oscillations 



that subsequently develop will occupy exactly the same space-time region for both the MNLS and defocusing 
NLS problems. A hint that this could be expected lies in the work of Kuvshinov and Lakhin [7] who 
found a mapping similar to ( 13 ) connecting the (genus one) Whitham modulation equations for MNLS and 
defocusing NLS. The latter are derived via an averaging procedure to model four coupled fields describing the 
bulk properties of a modulated periodic traveling wave solution, a structure that resembles the wave pattern 
initially generated by dispersive regularization of shocks. If one applies the formal matching procedure 
pioneered by Gurevich and Pitaevskii 0] to connect the Whitham description of a modulated wave with a 
background described by the dispersionless systems (T8T) and (10), one might expect to be able to formally 



establish that the explicit mapping ( 13 ) together with the mapping found by Kuvshinov and Lakhin would 



give rise to the same space-time region containing at least the first stage of development of the dispersive 
regularization process. The power of the observation that the same g-function controls both MNLS with 
supersonic initial data and also defocusing NLS in the semiclassical limit is that it makes the Gurevich- 
Pitaevskii procedure completely rigorous, and it also shows that similar results hold if shocks subsequently 
form in the solution of the Whitham modulation equations: no matter how many breaking curves or caustics 
appear in the semiclassical dynamics to separate regions modeled by multiphase oscillations of various genera, 
these caustics will always be exactly the same for both MNLS and defocusing NLS. 

The correspondence between the semiclassical MNLS equation with supersonic initial data and the defo- 
cusing NLS equation can be illustrated numerically. We solved the MNLS equation for a — 1 with initial 
data for which the Madelung variables were taken as: 



1 1 



-256x z 



and 



s(as,0) = l. 



(117) 



Note that since a = 1 we have Q e (x, 0) = a 2 p e (x,0) + au e (x, 0) — 1 = p e (x,0) > 1/10 so we indeed have 
a supersonic initial condition. We took e = 0.01, and we used a second-order Fourier split-step scheme 
[51 1151 [T7] that built in periodic boundary conditions on the domain x € [— 7r, it] with a spatial resolution of 
4096 Fourier harmonics and a time step of 5 • 10 -5 . Periodic boundary conditions are numerically convenient 
and an accurate representation of fixed boundary conditions if disturbances remain localized. We chose a 
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small positive initial background value of 1/10 for p because we wanted to be able to recover both p e and u e 
from the numerical data, and the latter requires careful interpretation if (f> e (x,t) becomes small. The split- 
step scheme alternately integrates the linear part of the equation exactly in the Fourier domain, and the 
nonlinear part by the midpoint rule in the physical domain using fixed-point iterations to ensure preservation 
of the L 2 -norm and hence maintain accuracy for integration over many time steps as is required to resolve 
the simulations when t< 1 (see [12 for a discussion of this important point). From p e (x,t) and u e (x,t) we 
calculated Q e (x,t) := a 2 p e (x,t) + au e (x,t) — 1 and then plotted p e (x,t) := p e (x,t)Q e (x,t). The resulting 
plot is shown in the left-hand panel of Figure [l] Now according to the correspondence established above 



MNLS tanh(20p) 



Defocusing NLS tanh(20p) 




Figure 1: Left: the function p e (x,t) recovered from a simulation of the MNLS equation with supersonic 



initial data consistent with (117), and with a = 1 and e = 0.01. Right: the function p e (x,t) recovered 



from a simulation of the defocusing NLS equation for e = 0.0122 with initial data obtained from (117) 



via the transformation (13) with a = 1. Dark color indicates small amplitude, while red (green) indicates 
positive (negative) values. Note that while as a square modulus p e (x, t) can never become negative under the 
defocusing NLS evolution, p e (x,t) can and indeed does become negative under MNLS evolution, although 
only after wave breaking occurs. Note that the oscillatory region appears to occupy the same spacetime 
region in both plots, exactly as predicted by the inverse-spectral theory. The Matlab code that generated 
these figures is available from the authors upon request. 

with the defocusing NLS equation in the semiclassical limit, we want to compare this simulation with one for 



the latter equation with initial data obtained in wavepacket form via the transformation ( 13 1 (with a = 1 
to match the MNLS simulation) of the Madelung initial conditions. We used virtually the same numerical 
method and chose e = 0.0122 (the slight discrepancy from the MNLS value of e ensures 27r-periodicity of 
the initial data for defocusing NLS). We computed p e (x,t) by taking the square modulus of the solution 
and our results are plotted in the right-hand panel of Figure [TJ One sees very clearly that the spacetime 
region in which the rapid oscillations appear (visible in both plots as a pattern of nearly parallel stripes) is 
almost the same for both plots. The theory outlined above indicates that the regions should agree exactly in 
the semiclassical limit, and also that outside of the (essentially common) oscillation region the plots should 
agree asymptotically. The green regions in the left-hand plot of Figure [T] are places where p t (x, t) has become 
negative, which is not inconsistent with our theory as long as it occurs only in the oscillation region for small 
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4 Subsonic Initial Data 



Now we assume that the initial data for the MNLS equation is globally subsonic for the dispersionless limit 
system: Q(x) := a 2 po(x) + auo(x) — 1 < for all x E K. Since a 2 po(x) > for all x and hence in this case 
1 — au(x) > for all x, it is obvious that the interval / defined in ( 85 ) is a subset of the positive real half-line 



Since in this case the discriminant of the quadratic on the right-hand side of (46) is negative for all 



x € R, the turning point curve is disjoint from the real axis in the complex z-plane. These facts imply that 
there are no real turning points x for any real value of z, and moreover, the reflection coefficient s(z, e) is 
small beyond all orders whenever z $ I, but is exponentially large for z € I. The "shadow" estimate proved 
in [2] cannot rule out a significant number of complex eigenvalues in this case, because the shadow of the 
turning point curve is an e-independent open subset of the complex z-plane. As pointed out earlier, finding 
both a leading-order estimate for the reflection coefficient in the interval / and good approximations for the 
discrete spectrum appear to require analyticity of (po,Uo) as functions of x. 

This situation should be compared with that arising in the study of semiclassical asymptotics for the 
initial- value problem for the focusing NLS equation with n = +1) with wavepacket initial data obtained 
from the subsonic data (po,i*o) via the mapping ( |14[ ). That is, we take initial data for focusing NLS in the 
form 



Mx) = M^y S{x), \ A(x) := y/- Po (x)Q(x), S(x) := S + / [u (y) + 2ap (y)] dy 



(118) 



where So is given by (88). To determine the asymptotic scattering data for the Zakharov-Shabat problem 



(27) with k = +1 corresponding to a wavepacket potential of this form, we make transformations directly 



analogous to ( 89 ) and ( 90 1 to arrive at the equivalent differential equation 



2ae- 



dW 



ZS.+l 

± 



dx 



iM zs ' +1 w| s ' + \ 



with coefficient matrix 



M 



zs. 



— 4z + 1 — au Q (x) — 2a 2 po(x) 
2ia v / -p {x)Q(x) 



-2iaty/-p Q (x)Q(x) 
4z — 1 + auo(x) + 2a 2 po(x) 



(119) 



(120) 



Seeking column vector solutions w = v(a;; z, e)e^ x ' z ^' / ( 2ae ) and applying the WKB method as in fjSJ one finds 
that M zs,+1 defined by (120) and M defined by (41) are similar matrices and so / satsifies again exactly 



the same characteristic equation (46) as in the MNLS case. The reflection coefficient 

s zs < +1 (z,e) :=r ZS ' +1 (A(z),e), zel 



(121) 



is therefore small beyond all orders for z ^ I but is large for z E I. If po and uo are analytic functions, then 
it may be possible (depending on the singularity structure of these functions in relation to the location of 
complex turning points) to express the asymptotic behavior of the reflection coefficient for z £ I in the form 



(z,e) - e ltE, W/ £ , where $ is given by the formula (|57[) in which X-(z) denotes an appropriate complex 



turning point (this makes 5R{i$(z)} > for z e /). In these cases the MNLS reflection coefficient s{z,e) is 



again indistinguishable from the focusing NLS reflection coefficient (121) in the semiclassical limit. 

Let us assume that / is an interval with a nonempty interior. There are indeed nontrivial cases where 
this is not so and the reflection coefficient is (essentially) everywhere negligible [5], but then one necessarily 
has 0(e~ 1 ) eigenvalues that contribute to the inverse-spectral Riemann-Hilbcrt problem, and the method we 
are about to explain to remove them from the problem does not apply. We suppose also that the scattering 
coefficients T\2{z,e) and 5 Z 2 S:+1 (A(z), e) extend analytically from I into a region R of the upper half-plane 
that contains all of the eigenvalues {zj} for the MNLS problem or {\a\j + 4} for the defocusing NLS 
problem. In this case, the proportionality constants corresponding to the eigenvalues are just given by 
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evaluating T\2{z 1 e) or S* 2 ^ ,+1 (A(z), e) at the eigenvalues. It follows that upon setting 



N^' +i (A(z)) 



N^ +i (A(z)) 



I s ZS, + l^ z e ^ e 2iS(z;x,t)/e 

1 



1 

_ S ZS, + 1( Z * e y e -2i0(z;x,t)/e 



N zs > +1 (A(z)), 



Z G R 



zeR* 



elsewhere, 



(122) 



and similarly 



NO) := < 



N(z) 



N(z) 



1 s(z,e)e 2ie(z ' x ' t ^ € 
1 



^s{z*,e)*e~ 2l0 ^ x ^/ e 



[N(z), 



z G i? 



z £ R* 



elsewhere, 



(123) 



the new matrix unknowns are both pole-free with jumps across the real axis (excluding positive parts of J) 
and the non-real boundary curves dR and dR* , both of which we take to be oriented from left-to-right, of 



R and R* respectively. Under suitable conditions, the approximate formula s(z, e) ~ s ,+ (z,e) ~ e 



i$(z)/e 

continues to hold for z G <9i? (this requires analytic continuation of the complex turning point x_ (z) into the 
complex z-plane) and we see that the asymptotic jump conditions for N(z) on dR and dR* take the form 



and 



N + (z)=N_(z) 



N+(z)=N_(z) 



I _ e i(20(z;x,t)+#O))/e 
1 



1 

_ e -i(2e(a;a!,*)+*(z))/e ^ 



Z G BR 



z G 0iZ* 



(124) 



(125) 



Here, the subscript "+" (respectively "— ") refers to the boundary value taken along the indicated oriented 
contour from the left (respectively right) side. Exactly the same asymptotic jump conditions hold on these 
curves with N(z) replaced by N zs,+1 (z). For the matrix N zs,+1 (z) there are only asymptotically negligible 
jumps across the real axis, while for the matrix N(z) this is true for z > while for z < one has the 
asymptotic jump condition N + (z) = z CT3 N_ (z)i~ CT3 . 

To study this problem in the semiclassical limit e <C 1, we can follow the methodology of Tovbis, Ve- 
nakides, and Zhou [16] by introducing a scalar function g(z), analytic for z G C\ (dRUdR*) and independent 
of e, and then setting, for the same ^-function, O(z) := N(z)e lg( - Z '>< 7 ^' i and O zs < +1 (z) := N zs ' +1 («)e i s , W <7!, / £ , 
which converts the jump conditions ( 124 1— ( 125 ) into the form 



and 



0+(*) = 0_(2) 



0+(z) = 0_(z) 



„i(g+{z)-g-(z))/e 


e i(g + (z)-g_(z))/e 
_ e -i(h+(z)+h-(z))/e 



s i{h + {z)+h-{z))/e 
-i{g + {z)-g-(z))/e 





-i{g + {z)-g-(z))/t 



z G dR 



z G dR*, 



(126) 



(127) 



where h(z) = 6(z; x,t) + h$(z) — g(z). Again, these also hold with O(z) replaced by O ZS ' +1 (z). The strategy 
used in [TH] is to choose g(z) = g(z*)* so that dR splits into a union of two types of subintervals: 
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• "Voids," subintervals in which g + — g_ = 0, while h+ + h— = 2h has a positive imaginary part. 

• "Bands," subintervals in which h + + is equal to a real constant, while g + — g_ is real and strictly 
increasing along dR in the direction of orientation. 

The jump matrix clearly converges rapidly to the identity matrix as e — > in the union of voids. In each 
band, a three-factor factorization and subsequent "steepest descent" lens deformation [TS] reduces the jump 
matrix to a constant (in z) off-diagonal form, with off-diagonal elements of unit modulus. 

We will not go further with the details here, but is should be clear that the use of the very same <?-function 
asymptotically reduces both the MNLS and focusing NLS Ricmann-Hilbert problems to similar limiting forms 
O(z) and O zs,+1 (z), respectively whose solutions can be obtained as in §3l with the only difference being 
that the band endpoints now are complex-valued and come in conjugate pairs. By estimating the errors 
carefully as in [TBI we may therefore simultaneously establish the validity for short times independent of e 



of the elliptic dispersionless limit system ( 10 ) for MNLS with subsonic initial data and also the validity of 



the elliptic dispersionless limit system QSJ) for focusing NLS with corresponding initial data provided via the 



mapping (14 1. As in the case of supersonic data for the MNLS equation, the subsonic MNLS and focusing 
NLS equations will display wave breaking phenomena in exactly the same regions of the (x,t) -plane in the 
semiclassical limit. 

Again we can illustrate these predictions with numerical simulations. Consider the initial data 

Pe(x,0) = ^q + I^ 256 * 2 and u € (x,0) = -l. (128) 

We take a — 1 and then it is easy to confirm that Q e (x,0) = a 2 p e (x,0) + au e (x,0) — 1 < — 7/5 < so 
this indeed corresponds to subsonic wavepacket initial data for the MNLS equation. The results of a Fourier 
split-step simulation on [— tt, it] (with 8192 Fourier modes and a time step of 3 • 10 -5 ) of the MNLS equation 



with initial data determined up to a constant phase by (128) are plotted in the variable p e (x,t) given by 
( 14 1 in the left-hand panel of Figure [2j Then, we solved the focusing NLS equation (using virtually the 



same method with 16384 Fourier modes and a time step of 10 ) with wavepacket initial data obtained from 



(1281 via the transformation (fMj), and a plot of the square modulus p e (x,t) of the solution is shown in the 



right-hand panel of Figure [2j The expected correspondence is indeed rather clear. 

5 Transsonic Initial Data 

As a specific example of transsonic initial data, consider 

Po(x) := sech 2 (x) and Uq(x) := 5 + /itanh(x), (129) 



A:= - „ and B := -f-. (130) 



for suitable real constants S and //. It is convenient to introduce in place of 5 and p. the following normalized 
parameters: 

. ? and B := — . 
4a 2 4a 

Let us assume that A and B are subject to the inequalities 

B 2 > A and B > \A\. (131) 

Then, it is easy to show that for such initial data, Q(x) — a 2 po(x) + au Q (x) — 1 is negative for x < x c and 
positive for x > x c with a simple root at x = x c = arctanh(2S — v 'AB 2 — AA + 1). Therefore, when t — 
we are in the subsonic case for x < x c and in the supersonic case for x > x c . 

In [3], the Cauchy initial- value problem for the MNLS equation subject to transsonic initial data of this 
type is carefully studied in the semiclassical limit e -C 1. It is easy to obtain specific information about the 
reflection coefficient s(z, e) from the WKB analysis presented in [J2] Indeed, only for z in the interval [zl, Zr] 
are there any real turning points at all, where 

z L :=a 2 (A-*-~^/l + 4B*-4A\ <z r :=0. (132) 
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Figure 2: Left: the function p e (x, t) recovered from a simulation of the MNLS equation with subsonic initial 
data consistent with (1281, and with a = 1 and e = 0.01. Right: the function p e (x,t) recovered from 



a simulation of the focusing NLS equation for e ~ 0.0078 with initial data obtained from ( 128 ) via the 



transformation (14) with a — 1. Dark color indicates small amplitude, while red (green) indicates positive 
(negative) values. As in Figure [l] we see that p t {x,t) can become negative under the MNLS evolution, and 
also that the oscillations approximately occupy the same spacetime region in both cases. The Matlab code 
that generated these figures is available from the authors upon request. 



Outside of this interval, s(z, e) is either exponentially small or exponentially large. Since m± := 5 ± ji in 
the present situation, according to the discussion around the definition (85) of the interval /, the latter case 
only occurs for < z < a 2 (A + B). 

For both z\, < z < and < z < a 2 (A + B) the reflection coefficient is given asymptotically by the 
formula s(z, e) = e l *^' e up to a relative error that is of order O(e). In [3] this approximation is made 
completely rigorous in a direct fashion; rather than using the WKB method, one can (only for initial data of 
the special form (129)) solve the direct spectral problem (15) exactly in terms of hypergeometric functions 
[2J. This allows the expression of all scattering data for all e > in terms of gamma functions without 
approximation, which in turn allows the semiclassical approximation of scattering data to be carried out in 
full detail with the use of Stirling's formula and its variants. One of the outcomes of this approach is that 
there are no eigenvalues whatsoever for any e > as a consequence of the first of the two inequalities ( 131 ). 
The function $(2) defined generally by the WKB formula (57) turns out to have an analytic continuation 
from the interval (zl, a 2 (A + B)) to the whole upper half z-plane. As a real- valued function in the subinterval 
(zl,zr), this analytic function in turn can be continued to the domain z G C \ ((— 00, zi] U [zr, +00)) as 
a function satisfying the Schwartz-symmetry condition $(z*) = $(z)*. The latter function has the explicit 
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expression (obtained either from (57 1 or via Stirling asymptotics): 



$0) = -4aSlog(2) - — (a 2 (A + B) — z) log i-(a 2 (A + B) - z) 
a \a 



2 ({a 2 B 2 - z) 1 ' 2 + aB) log (2 [{a 2 B 2 



z 



1 / 2 +aB 



+ 2(a 2 B 2 - z) 



aB)) 




z) 1 ' 2 - 


a 


z) 1 ' 2 + 


% 

a 



(133) 



Here all logarithms and square roots refer to the principal branches. The exponent of the reflection coefficient 
is obtained from this formula by direct evaluation for z G (zl, zr = 0) and by taking a boundary value from 
the upper half-plane for z € (zr, = 0, a 2 (A + B)). The latter boundary value satisfies: 

'2tt {{a 2 B 2 - z) 1/2 - aB) , < z < z+ 

M* <->} = < (134) 

— (z-a 2 (A + B)). z+ < z <a 2 (A + B). 

a 



Here 



a 2 ( A- i + ^y/l + AB 2 -AA ) e (0,a 2 (A + B)), (135) 



and we note that indeed 3{<i>+(2;)} < for < z < a 2 (A + B) as is consistent with s{z, e) ~ e l ®+( z )/ f - being 
exponentially large. 

The semiclassical asymptotic analysis of the Riemann-Hilbcrt problem for MNLS described in §1.2.1 



necessarily combines simultaneously essential elements of both the type of methodology described in <|3_ 
(the case of globally supersonic initial data) and also the type of methodology described in ^j4] (the case of 
globally subsonic initial data). One has to factor the jump matrix into two factors in an interval of the 
form (w, a 2 (A + B)) with Z\, < w < zji = 0, an interval containing that in which the reflection coefficient is 



large, and deform the two factors into the upper and lower half-planes by a transformation like (123). Here 
one has to be careful because while there are no eigenvalues, there exist poles of the scattering coefficient 
Txziz, e) emanating from the point z = a 2 (A + B) into the upper half-plane, and the transformation (123) 



will introduce artificial poles (in [5] they are called "phantom poles" , but they are also similar to "false poles" 
[13] or "resonance poles" [14] in scattering theory). A similar situation occurs in the semiclassical analysis 
of focusing NLS in [16], and the authors of that work argue that the region R can be chosen to abut the 
line of poles only at the right endpoint of the interval in which the reflection coefficient is large; this reduces 
the problem to a local one near the endpoint, and there an implicit parametrix is constructed to deal with 
the breakdown of the semiclassical approximation of the reflection coefficient. By contrast, in [5] a local 
construction is used that actually allows poles to reside in R (by a modified pole-removing substitution) and 
controls all errors with exponential accuracy near z — a 2 (A + B) without the use of any local parametrix 
whatsoever. 

Once the large reflection coefficient has been removed from the real axis in this way, one introduces a 
g-function that now has the job of controlling at the same time the triangular jump matrix factors on dR 
and OR* as in ^4] and the oscillatory full jump matrix in the real interval (zl,2r = 0) as in Sj3] One of 
the results of [3] is that the transsonic point x c actually persists for small nonzero t as a moving boundary 
x = x c (t) that separates asymptotically subsonic flow for x < x c (t) from asymptotically supersonic flow for 
x > x c (t). In the gas dynamics literature, such a moving boundary is called a sonic line. At the level of the 
Riemann-Hilbert problem and its g-function, what happens as one crosses the sonic line from left-to-right is 
that the complex endpoint q = q(x,t) of a single band on dR connecting z = q to z = a 2 (A + B) collides 
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with its complex conjugate at a point z c (t) £ (zl, zr) and the pair re-opens as a real band (ao(x, t), bo(x, t)) 
surrounded by voids; the point w at which dR meets the real axis in the interval (zl, zr) is then taken as the 
right cndpoint z = &o of the newly-born real band. To determine the sonic line, one first solves the equation 

m *>»( ^ 2 r 2(A+B) Wtoldf n 

__$ (Zc) + _y o __ =0 (136 ) 



for z c = z c (t) (for i = this gives z c = a 2 (2B 2 — A — y '4B 2 — 4^4 + 1) e (zl,0), which can be used as 
a starting point for Newton iteration to obtain the solution for small nonzero t), and then the sonic line 
x = x c (t) is given explicitly by 

* W -id- 4„(„, . + f - f /" V+fl) =*^. (137) 

a 4 27r J s-« c (i) 

In the vicinity of the sonic line, and in fact throughout the subsonic flow region to its left, the e-dependent 
functions p e (x, t) and u e (x, t) are accurately approximated in the limit e -C 1 by a certain solution p(x, t) and 
u{x,t) of the mixed- type dispersionless MNLS system (10). Moreover, the solution is characterized without 
dealing with partial differential equations at all by the following formulae. For arbitrary p > and 
let S(z; p, u) be the function determined by the conditions that 

165(z; p, u) 2 = (4z - 1 + au) 2 + 16a 2 pz (138) 

and that S{z; p, u) is analytic in the complex plane excluding a branch cut along the straight line connecting 
the two roots of the quadratic ( 138 ), with asymptotic behavior S(z; p,u) = z + 0(1) as z -> 00. Then define 

. 27T, , 1 _[ ®( S )ds f a2(A+B ^{^'{ S )}ds 

M {p,u;x,t) := -^{ax- (au + 2a p)t) - — f — r + / — — 

a 2 4i / S(s;p,u) J S{s;p,u) 

Mi (p, u; x,t) := — ^ ((1 — au — 2a 2 p)ax + (6a 4 p 2 + 6a 3 pu — 4a 2 p + a 2 u 2 — au) t) (139) 

1 r<5>'{s)sds r a2 ( A + B ^{& + { s )}sd s 



Ai J S(s;p,u) J S(s;p,u) 

Here, the contour integrals are taken over a positively-oriented loop surrounding the branch cut of S. One 
then solves the equations Mo = Mi — by Newton iteration for (p, u) as functions of (x, t) near t — with 
the initial conditions p(x, 0) = po(x) and u(x, 0) = u (x). The Jacobian of this system of equations vanishes 
along the sonic line, even though the solution p{x,t) and u(x,t) remains smooth. This means that Newton 
iteration to find (fi,u) will fail for x — x c (t), and will be numerically difficult (will require a very accurate 
initial guess for convergence) for (x, t) close to the sonic line x = x c (t). However, along the sonic line exactly, 
p(x c (t),t) and u(x c (t),t) can be found directly from the auxiliary function z c (t) solving (136): 

p(x (t),t) = —=Zc(t) and u(x c (t),t) = -(l + 4z c (t)). (140) 
a z a 

It is easy to check that upon elimination of z c between these two equations one arrives at the condition 
Q = a 2 p + au — 1 = 0, the condition for degeneration of the type of the quasilinear system ( 10 ). 
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